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A perturbative method for solving the Langevin equation of inflationary cosmology in presence 
of backreaction is presented. In the Gaussian approximation, the method permits an explicit calcu- 
lation of the probability distribution of the inflaton field for an arbitrary potential, with or without 
the volume effects taken into account. The perturbative method is then applied to various concrete 
models namely large field, small field, hybrid and running mass inflation. New results on the stochas- 
tic behavior of the inflaton field in those models are obtained. In particular, it is confirmed that 
the stochastic effects can be important in new inflation while it is demonstrated they are negligible 
in (vacuum dominated) hybrid inflation. The case of stochastic running mass inflation is discussed 
in some details and it is argued that quantum effects blur the distinction between the four classical 
versions of this model. It is also shown that the self-reproducing regime is likely to be important in 
this case. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

Quantum effects play a crucial role during inflation. In 
particular, they are responsible for the self-reproducing 
behavior of the Universe also known as "eternal infla- 
tion" 0. In this regime, the quantum fluctuations are 
so important that they can dominate the classical dy- 
namics and, therefore, treating them properly becomes 
mandatory. Technically, this is a difficult task especially 
when it is necessary to take into account the backreac- 
tion of the quantum field on the geometry. Stochastic 
inflation 0, Q, 0> S IS IE aims at providing a for- 
malism where the previous difficulties can be partially 
circumvented. In the stochastic inflation approach, one 
is mainly interested in the evolution of a coarse-grained 
field, typically the original scalar field averaged over a 
Hubble patch, and the quantum effects are modeled by a 
stochastic noise originating from the small-scale Fourier 
modes. Consequently, the dynamics of the coarse-grained 
field is controlled by a Langevin equation. Then, en- 
dowed with a solution of this equation, one can compute 
the probability density function of the field and the var- 
ious correlation functions. 

Even if stochastic inflation simplifies the calculation 
of the quantum effects, the Langevin equation remains 
difficult to solve (without, of course, relying on numerical 
computations). The case without backreaction (in a de 
Sitter background) has been investigated in Ref. J8] where 
it has been shown that solutions for an arbitrary potential 
can be obtained. The case with backreaction is clearly 
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much more complicated and something can be said about 
the solution only for very specific potentials. The usual 
approach is applicable when the inflaton potential is such 
that an exact solution exists for some power of the field. 
Then, in order to obtain the stochastic field itself and 
its various correlation functions, an expansion in terms 
of the coupling constant of the potential is performed. 
Explicit examples of this approach can be, for instance, 
found in Refs. 0, and more recently in Ref. an d 
are briefly discussed in the following. The main point 
is that it is first necessary to obtain a solution to be 
subsequently able to perform the expansion. This is the 
reason why the applicability of the method is severely 
limited. 

In this article, we present a method based on a pertur- 
bative expansion in the stochastic noise, the expansion 
being performed directly in the Langevin equation. As 
a consequence, our method does not require getting first 
a solution and, a priori, can be pushed to any order, 
the only limitation being the mathematical complexity 
of the obtained expressions. This represents a crucial 
advantage over the other approaches which allows us to 
treat analytically the case of an arbitrary potential with 
backreaction. It should be noticed that the idea to solve 
perturbatively the Langevin equation was put forward for 
the first time in Ref. [IJ. In that reference, the method 
was used to compute the three-point correlation function 
of the Cosmic Microwave Background (CMB) anisotropy 
while, here, we use it in order to study the influence of the 
quantum effects on the behavior of the background infla- 
ton. In particular, we show that, at second order, the 
calculation of the probability density function reduces 
to the calculation of a single quadrature. Moreover, the 
case where the volume effects are taken into account only 
requires the calculation of an additional integral. The 
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method is then applied to various concrete cases, as the 
chaotic, new, hybrid and running mass inflationary mod- 
els. In each case, the probability density function can 
be computed analytically and the volume effects evalu- 
ated exactly. This allows us to study the relevance of 
the quantum effects in those models. In particular, in 
the case of the running mass model, this is the first time 
that such an investigation is carried out. 

Our paper is organized as follows. In the next section 
(Sec.[nj, we briefly review the basic equations of stochas- 
tic inflation. Then, in Sec. II I II we present our method 
and compare it with the approaches known in the existing 
literature. We also show how to compute the probability 
density function, with or without the volume effects taken 
into account, directly from the Langevin equation with- 
out writing a Fokker-Planck equation. As already men- 
tioned, we explicitly demonstrate that this calculation 
simply reduces to the calculation of a single quadrature. 
In Sec. lIVI we briefly present the inflationary models the 
stochastic effects of which are computed in the subse- 
quent section. In particular, we focus on the choice of 
the free parameters characterizing the corresponding po- 
tentials since their values are crucial in order to estimate 
the importance of the quantum effects. In Sec.[V] we dis- 
cuss and interpret our results for the cases of large field 
models, small field models, hybrid inflation and running 
mass inflation. To our knowledge, the calculation of the 
mean value and the variance of the coarse-grained field 
in the last three models was never done before (when the 
backreaction is taken into account). We end this article 
by some concluding remarks. 



II. BASIC EQUATIONS 

In the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) Universe, the assumptions of homogeneity 
and isotropy allow us to write the metric in the simple 
form ds 2 = — dt 2 + a 2 [t)8ijdx i dx^ , where a(t) is the 
time-dependent scale factor and where we have assumed 
flat space-like sections. In such a space-time the evo- 
lution of an homogeneous scalar field <j>(t), sourcing 
the metric evolution, is described by the Klein-Gordon 
equation 



4> + 3H<j> + V'{(f>) =0, 



(1) 



where a dot means a derivation with respect to the cos- 
mic time t and a prime a derivation with respect to the 
scalar field <\>. This equation is coupled to the Friedmann 
equation for the scale factor 



H z = [ - 

a 



(2) 



where we have defined n = 87r/r 



being the Planck 



mass. 

During a phase of inflation, if the slow roll approxi- 
mation is satisfied then the acceleration 6 of the field is 



negligible compared to the friction term 3H<j) and, at the 
same time, the kinetic energy 4> 2 /2 is small compared to 
the potential energy V((j>). This approximation consider- 
ably simplifies the equations describing the evolution of 
the system. These ones can now be re-written as 



ff 2 W>) = -i/W>), 4> + -h'{^) =o. 



(3) 



At the classical level, there is nothing more to say. Once 
we are given a potential V(4>), the above equations can 
be solved and the time evolution of the scale factor and 
of the infiaton field obtained. 

The problem becomes much more complicated when 
the field <f> is considered as a quantum operator. A first 
difficulty arises because quantizing a scalar in curved 
space-time is technically complicated for the case of an 
arbitrary potential (for an arbitrary potential, the Klein- 
Gordon equation is non-linear and one cannot Fourier ex- 
pand the field and write an equation for a time-dependent 
mode function). A second (more fundamental) diffi- 
culty is to take into account the effect of the quantum 
scalar field on the geometry, i.e. the backreaction prob- 
lem. Since the infiaton sources the Einstein equations, 
the Friedmann equation (taken literally) indicates that 
the geometry should also be quantized. Unfortunately, 
this quantum gravity regime is presently not known and 
the previous program cannot be carried out. 

The stochastic formalism allows us to circumvent these 
difficulties. In the stochastic formalism 0, one is inter- 
ested in the dynamics of a "coarse-grained" field ip(t, x). 
This coarse-grained field is defined to be the spatial av- 
erage of the ordinary field <fi over a physical volume the 
size of which is typically larger than the Hubble radius 
H^ 1 = a/a. Therefore, ip(t,x) basically contains the 
long- wavelength Fourier modes (i.e. those with comov- 
ing wavenumber such that k < aH) only. 

The evolution of the coarse-grained field is still de- 
scribed by the Klein-Gordon equation but a suitable 
random noise field acting as a classical stochastic 
source term, should be added in order to take into ac- 
count the effect of the quantum fluctuations. In the slow- 
roll approximation, the evolution of the coarse-grained 
field is thus governed by a first order LangOTin^lfice^iffer_- 
ential equation which can be written as [lj, 



1 dV H 3 / 2 „ . 



(4) 



where the noise field £ is defined in such a way that its 
mean and two-point correlation function simply read 



<£(*)> = 0, mm)=s(t-t'), 



(5) 



8(z) being the Dirac distribution. The normalization of 
the correlation function is chosen in order to reproduce 
the ordinary result, (y 2 ) = H 3 t/(4:ir 2 ), valid for a free 
field in de Sitter space-time. 
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It is at this point that the backreaction problem shows 
up. The standard assumption is that the Hubble param- 
eter in Eq. is only controlled by the coarse-grained 
field. Then, one needs to specify how H depends on 
ip and one naturally assumes that the Friedmann equa- 
tion (in the slow-roll approximation) holds for the coarse 
grained quantities, namely 



H\ip)~-V{ V ). 



(6) 



A direct consequence of the above equation is that the 
noise becomes multiplicative. The formalism briefly de- 
scribed previously also indicates that the coarse-grained 
field <p now describes a Brownian motion for which the 
classical drift is modified by the quantum diffusion term. 
The main goal is to solve Eq. since, endowed with 
the solution, we can then evaluate the probability den- 
sity function of the coarse-grained field and/or various 
correlation functions. 



III. SOLVING STOCHASTIC INFLATION 

A. Perturbative Method 

As discussed above, the main purpose of this article is 
to present and study a method for solving the Langevin 
equation. This method was used for the first time in 
Ref. ^3 f° r the calculation of CMB non-Gaussianities 
and in Ref. |13J in order to compute how the quantum 
effects affect the behavior of the quintessence field during 
inflation. Here, we develop the method in full generality 
including the calculation of the probability density func- 
tion. The main idea is to consider the coarse-grained field 
ip as a perturbation of the classical solution ip c \ (the clas- 
sical solution is defined as the solution of the Langevin 
equation without the noise), (p c \ being supposed to be 
known. The corrections to tp c \ are obtained by adding 
successive terms of higher and higher powers in the noise, 
i.e. 



(7) 



where the term 8<pi (t) depends on the noise at the power 
i. The equations of motion controlling the evolution of 
the term 6<fi(t) are obtained by inserting the above ex- 
pansion into the Langevin equation and by identifying 
the terms of same order. Expanding up to second order, 
we get two linear differential equations for Sipi and 5<f2, 
namely 



— — + -H {ipci)5(pi = f (t) . 

dr K 2tt 



(8) 



and 



---- - 7// "^Vf 



— — + -H {ip c i)5ip2 = - 
at k k 



(9) 



These equations are similar to Eq. (17) of Ref. [12(- The 
only difference is that, in the previous reference, the 
Langevin equation is written in term of the number of 
e-folds while, here, our time variable is the cosmic time. 
Since the above equations are linear, they can be solved 
by varying the integration constant. Straightforward ma- 
nipulations lead to 

%(1) ^f jT £!m H , (10) 



2tt 



t iD H> [ 9c1 (t)} 



where we have assumed that the initial conditions are 
such that 6<pi(t = t- ln ) = 0. In the same manner, the 
solution for 8<p 2 (t) can be easily obtained and reads 

TTl ft TT 11/ 

SMt)=-~ Jf dr— 8 V \{r) 

+ ^£dT# 1 / 2 <Vi(rX(T). (11) 

As expected, 8<pi is linear in the noise £ while 8ip 2 is 
quadratic. Of course, the expansion could be pushed 
further and one could evaluate Sips, 8tp4 etc . . . using the 
same technique. 

We are now in a position where the various correlation 
functions can be calculated exactly. Since 8(p± is linear 
in the noise, its mean value obviously vanishes 



(8(fi! 



0. 



(12) 



This means that Sipi does not introduce any correction to 
the mean value of the coarse-grained field. As a matter 
of fact, Sipi directly contributes only to the variance of 
(p. Using the white noise correlation function given by 
Eq. JSJ, we obtain 



2 V2^ 



Vol 



(13) 



We see that the calculation of the variance reduces to a 
simple quadrature. In order to calculate the correction 
to the mean value we must consider the mean of 6if2- 
Using the fact that (8<Pi(t)£(t)) = H 3 / 2 /(An), we arrive 
at 



(Sip 2 





Jtp ci 


2irm 2 
pi 










H" ( 


•>Vcl 







H_ 
~H~< 



3/P 
~2W 



(14) 



In this expression, the first term is nothing but the one 
given in (JT3J, i.e. (8<Pi), while the second one can be 
evaluated exactly with the help of an integration by parts. 
This leads to 



(8<P2 



IT 
2H 



7 



H' 

Airm 2 



Hi 



H 3 

m 2 



(15) 
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Therefore, at second order in the noise, everything can 
be reduced to the calculation of a single quadrature, the 
one of Eq. (fH5|) . Before discussing the probability density 
function, we compare the method described above with 
what is already known in the literature. 



B. Comparison with Other Methods 

As already mentioned before, in order to treat the case 
with backreaction, various papers 0, IToj concentrate on 
very particular cases where the Langevin equation can 
be solved exactly (see also the recent article, Ref. [Tlj . 
where the same method is used) . The typical example of 
this procedure is the model described by a quartic poten- 
tial, V = 3A4((5 4 /(87r), A4 being a dimensionless coupling 
constant, for which the Langevin equation can be written 
as 



d 
dt 



2tt 



■if : 



2vr 



Jf_ 

m D1 



m- (is) 



This equation can be solved exactly because it takes the 
form of a Bernoulli equation after a change of variable. 
However, in this case, one does not obtain the coarse- 
grained field itself but rather some power of it, namely 



¥T 2 (f)=^ 2 (*) [!-¥(*)] , 



(17) 



where ip c \(t) is the classical solution (which, in the slow- 
roll approximation, is known explicitly) and where the 
stochastic quantity ^(t) is defined by 



*(*) 



\3/4 1/2 
A 4 "Vi 



dr 



1 2 



£(r), (18) 



which is a new dimensionless Gaussian noise with van- 
ishing mean value and whose variance (and higher cor- 
relation functions) can easily be computed. Therefore, if 
one wants to obtain the field itself, it is necessary to take 
the inverse square root of the solution l|17(l . 

At this point, several remarks are in order. First, the 
only way to compute the coarse-grained field ip and its 
various correlation functions is to expand (1 — vf) -1 / 2 
in 'J, that is to say in the coupling constant A4, and to 
truncate the expansion at some order (the series does 
not converge anyway). Thus, we see that, despite the 
fact that we have an exact solution, an expansion is still 
required in order to use Eq. lfT7|) concretely. Second, one 
can show that the expansion in the coupling constant is 
equivalent to our expansion in the noise |l3j . However, 
clearly, our method is more general because it is not re- 
stricted to the situation where an exact solution of the 
Langevin equation is available. This is because the ex- 
pansion is directly performed in the Langevin equation 
rather than in its solution. The drawback of the method 
used in Refs. is clearly that it is first necessary to 

find a solution of the Langevin equation before the ex- 
pansion can be taken. We will illustrate this last remark 



on the calculation of the criterion which determines when 
the self-reproducing regime becomes efficient (i.e. when 
the quantum fluctuations dominate the classical drift). 
In Ref. ^l|, using the model V (x ip 4 for the reasons 
described before, the authors have recovered the stan- 
dard result that this happens when the initial value of 
the inflaton is larger than ip in ~ A 4 m pl . Using our 
formalism, we will derive this criterion for any model of 
the form V oc tp n , which would not have been possible 
with the method of Ref. . 

Another possibility studied in the literature is the so- 
called scaling solutions method, see Ref. The idea 
is to perform a change of variable and to work in terms 
of a new stochastic process r](tp) such that the Langevin 
equation takes the form drj/dt = J-(t; rj)^(t), where T is 
a priori a complicated function of 77. If rj is replaced by 
?/ c i in J 7 , then the new Langevin equation becomes solv- 
able. Therefore, one sees that this method bears some 
resemblance with the method investigated here. How- 
ever, there also exists important differences. First, our 
method does not require any change of variable which is 
an advantage since, in general, the link between ip and rj is 
quite complicated. Second, in the scaling method, there 
is no systematic expansion in the noise (in some sense one 
always works at first order) while in our method we can 
go to any order, the only limitation being mathematical 
complexity. Third, a saddle point approximation is used 
to estimate the effective dispersion while, in our case, 
following the calculations of the previous subsection, this 
can be done exactly. Fourth, it is difficult to evaluate the 
reliability of the scaling limit while we will discuss in a 
forthcoming article [l5l ] how the accuracy of our method 
can be determined precisely. Finally, let us also stress 
that, in Ref. Q, only the cases of large field models 
and exponential potentials are considered (in principle, 
it would be possible to treat other models with the scaling 
method although it is unclear whether this would lead to 
analytical expressions for the probability density) while, 
here, we will also apply our method to the new, hybrid 
and running mass inflationary scenarios. 

Finally, let us repeat that the method studied here is 
similar to the one used in Ref. even if we apply it in 
a different context. In particular, the solution given by 
Eqs. itTUll and HTJ are identical to Eqs. (19) and (21) of 
Ref. [I2j, these formulas, however, being written in terms 
of the total number of e-folds rather than in terms of 
cosmic time. In Ref. |l^. these results are applied to the 
calculation of the three-point correlation function of the 
CMB fluctuations while, in the present article, we use 
them, among others, in order to derive the probability 
density function of the background field. 

To end this subsection, let us recall that more details 
on the other methods discussed here can be found in 
Ref. [H. 
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C. Probability Density Function 

After having compared our approach with other for- 
malisms, we now come back to the general method and 
show how one can calculate the single point probabil- 
ity distribution P c (ip,t) of the coarse-grained field (also 
sometimes called P p in the literature) . Let us recall that 
P c (ip,t) is the probability of the stochastic process to 
assume a given value at a given time in a single coarse- 
grained domain. Very often, this probability distribution 
is obtained from a Fokkcr-Planck equation. However, it 
can also be determined from U&, 



p (<p,t) = 



(19) 



with the normalization Ao simply being given by A/"o = 
(JV^e-^ ' 1 ^ 2 )- 1 and where we have introduced the 
definition f T g = J dr/(r)<7(r). As one can see on the 
above expression, the noise probability distribution is 
Gaussian and correctly yields the noise correlation func- 
tion (£(*)£(*')) = C(t,t'). In Eq. (HHJ), the stochastic 
field will be given by our perturbative solution, namely 
</?[£] = ip c \ + Sipi + 6if2- The first and second order cor- 
rections to the classical solution, which are linear and 
bilinear in the noise respectively, can also be written as 
5<pi = J T £ and Sf2 — ? T A£, where the detailed defini- 
tion of J and A are given in Appendix lAl Then, using 
the integral representation of the S function, we get 



where tp[£] is the solution of the Langevin equation, and 
the mean value has to be evaluated with the functional 
probability distribution V[(] of the noise 



V[i] = A" exp 



1 



(20) 



P c ((fi,t) = ^- / dyexp[iy(ip cl - ip))J\f / £>£ exp 



(21) 



In this expression the functional integration is a Gaussian integration involving the kernel C 1 — 2iyA. Defining the 
new coefficient Af y ~ e~^ ( c ~ 2l J' A )C/ 2 ]~ 1 j the functional integration yields 



P C (V, t) = ^ J dy exp [iy (<p cl - <p)\ j^- exp 



-~y 2 J T (C- 1 -2iyA)- 1 J 



(22) 



r 



Then, generalizing the relation valid for two finite n x n 
matrices M and N, 



< \n xe -x i M iJ Xj 



d".rc '- v r 



detN 
detM 



_ e -TrlnN- 1 M/2 



(23) 



to the continuous case 01 1 we can wr he the ratio of the 
two normalization coefficients as 



Ao 



exp 



--Trln(l - 2iyAC) 



(24) 



up to second order in the noise. Finally, if the argument 
of the second exponential in Eq. (1221 is expanded in pow- 
ers of the noise, then all the terms but J T C J = (Sipf) can 
be neglected. Evaluating the remaining ordinary integra- 
tion over y we get the normalized Gaussian distribution 



1 



: exp 



2 1 



2(S<pl) 



(25) 



This distribution is centered over the mean value (tp) ~ 
(f c \ + (5ip2) with variance \8ip\j. The above equation is 



one of the main result of this article. In order to evaluate 
P c only the integration (|13fl is necessary which illustrates 
the power of the perturbative method. 



D. Volume Effects 

If we want to investigate the evolution of the probabil- 
ity distribution of the field when spatially averaged over 
the entire Universe (and not only in the single domain) 
we must take into account the volume effects. These are 
due to the fact that the size of each homogeneous domain 
depends on the value of the field within the domain itself, 
and we expect larger domains to give a more important 
contribution to the average over space. The value of each 
field-dependent quantity must thus be weighted with the 
physical volume a 3 (ip) = exp[3 J drH((p)]. Therefore, 
the normalized probability distribution accounting for 
volume effects can be obtained from 



Py(<P,t) 



(e 3 / 



(26) 
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In order to compute this new distribution function, we 
expand perturbatively ff (<£>[£]) up to second order in the 
noise, 



I T £ while the second order 
j!dr [H>(t)6Mt)+H"(t)5<pI/2 



ones take 

Tt 



the form 



H(<p[£]) = H cl + H'JSfi + S<P*) + ^fW? , (27) 

and write the first order term appearing in the argument 
of the exponential in Eq. il26|) as J t drH' (t)5(P\(t) = 



Then, one can repeat the calculations performed in the 
previous subsection and we obtain the following expres- 
sion for the volume weighted distribution function 



Pv(<P,t) 



2ir (a 3 (ip)} 



dy exp [iy (<p cl - ip)] M / T>£ exp 



-^(C- 1 -2iyA-6B)t+(iyJ 1 



3/ T )? 



(28) 



This is again a Gaussian integration but with the modi- 
fied kernel C^ 1 — 2iyA— 6B where the term — 6B is a new 
contribution that accounts for the volume effects. These 
volume effects also manifest themselves in the term linear 
in the noise (i.e. the term proportional to 7). Following 



the same steps as before, we define the new normaliza- 

-l 

and evaluate 



tion My = 



the functional integral. We obtain 



J 



Py(<P, t) 



1 



2tt (a 3 (ip)) 



dy exp [iy (ip c 



NT Mo 

<P)] Jp exp 



~(iyJ T + 3/ T ) (CT 1 - 2iyA - 6B) 1 (iyj + 31) 



(29) 



As in the previous case, the functional inverse of the mod- 
ified kernel (C _1 — 2iyA — 6B) _1 simply reduces, up to 
second order, to C. 

We must now evaluate the denominator. Exactly in 
the same way as before, this term becomes 

(a 3 (p))=M Jv^e^ V -i ? T (c -i_ 6B)e + 3/ T c ) (30) 

and the functional integration yields (^(ip)') ~ 
(Mq/Mq) exp (9/ T C//2), where Mq is a new normaliza- 
tion coefficient the explicit expression of which we do not 
give here for simplicity. Then, we insert this last expres- 
sion in Eq. J22J|. As expected the terms Mo cancels out. 
Moreover, when evaluating the ratio Mq/ML the contri- 
butions of the terms describing the volume effects (those 
involving B) also cancel out at second order and we get 
No/My ~ (J-y^vi) as before. Putting everything together, 
we obtain the final result 



exp 



(p-(p)-3lVy 
2{6&) 



(31) 



where (tp) = ip c \ + (5tpi) is the usual mean value. 

This expression should be compared with Eq. 1)25(1. 
The variance of the resulting Gaussian probability dis- 
tribution is thus unchanged, while the volume-weighted 



mean value (ip) = (p) + 3 I T J gets the extra correction 



3/ T J = 3 



(32) 



'in 



So far, the calculation of the volume effects, in partic- 
ular Eqs. 113111 and l|32|) . do not rely on the slow-roll ap- 
proximation. However, if this approximation is satisfied, 
then the volume contribution can be easily calculated for 
a generic potential. We obtain 



3/V = 



12ff '^d^-12^^M. (33) 



Therefore, we see that the calculation of the volume ef- 
fects only requires the computation of one additional 
quadrature. 



IV. INFLATIONARY MODELS 

In this section, we briefly present the inflationary mod- 
els to which our method is applied in the next section. 
In particular, we carefully discuss the choice of the free 
parameters characterizing those models since their nu- 
merical values turn out to be crucial in order to estimate 
the importance of the stochastic effects. 



A. Generalities 

We adopt a parameterization of V(<p) suitable for de- 
scribing different types of inflationary models, and we 
write the potential as 



V((p) = M 4 



(34) 



where a = 0, I and b = ±1 according to the case under 
consideration, while M, y, and n (with n > 2) are free 
parameters. If a = and b — 1 we have monomial po- 
tentials describing chaotic inflation (T^i also commonly 
known as "large field models" (LF) because the initial 
value of the field (rolling towards the origin) is typi- 
cally much larger than m pj . The case n — 4 has al- 
ready been treated in Refs. 0,0] while the general case, 
i.e. for an arbitrary value of n, was studied for instance 
in Refs. 0,0, El- Quantum effects for potentials with 
a = 1 have not been computed explicitly before and, 
therefore, we will mainly focus on those examples. Po- 
tentials with a — 1 and b — — 1 belongs to the class of 
the "small fields models" (SF) such as the new inflation 
scenario, where the field starts in the false vacuum close 
to the origin and moves down to <fi = fi as in a sponta- 
neous symmetry breaking |19| . At this point, one remark 
is in order. In fact, the case of stochastic new inflation 
has been investigated many times in the literature, for 
instance in Ref. 4]. But usually, and this is in this sense 
that the treatment presented here is new, the backreac- 
tion is not taken into account and the Hubble parameter 
is just considered as a constant. In this article, we do not 
make this assumption. Finally, the case a = 1 and b = 1 
describes hybrid inflation pol ] . Although hybrid inflation 
is a two-field model, the slow-rolling phase taking place 
in the inflationary valley of the potential can effectively 
be described as a single field model. 

We also consider the running-mass model (RM) 0, 
l22l | , the potential of which does not belong to the class 
presented above. For this model V is given by [2l| 



V{<p) = M 4 



In 



Ml 



2 1 



(35) 



where M pl = m pi /\/8Tr. In this expression, M, c and tpo 
are free parameters (In Ref. |2lj . M 4 is denoted Vq and 
(po is written 0* ) . Let us notice that c can be positive or 
negative. 

Our next step consists in obtaining the classical trajec- 
tory for these models. This can be done if the slow-roll 
approximation is satisfied but, even in this case, the clas- 
sical trajectory can be found implicitly only. In terms of 



total number of e-folds N, we have for the models de- 
scribed by Eq. 



N 



2 f<Pcl/V 

-/A I dxx 1 -" {a + bx n ) . (36) 



no 



The integration can easily be performed and the solution 
can be expressed as (in the following we use the fact that, 
when non-vanishing, a is one and that b is just a sign) 



N = K 




2-n 



(37) 



for n 7^ 2, while for n = 2 one has 



N = K- 



ab\n 



<gcl 



(38) 



If a = it is very easy to find the field evolution inverting 
the above expressions and solving for 4>. This was done, 
for instance, in Ref. 01 . On the contrary, if a = 1 an 
explicit solution can be found only for particular values 
of n. For simplicity, we concentrate on the specific case 
71 = 2 for which we get 



Yd 



bW 



b 



exp 



<pl - (N/2n)ml 
bfi 2 



(39) 



where Wq(x) is the principal branch of the Lambert func- 
tion (2^. This special function is the solution of the 
equation W(x)e w ^ = x. Since the curve xe x has a 
global minimum for x = —I, its inverse is a multivalued 
function with two branches on the real axe (and infinite 
branches on the complex plane). The one being contin- 
uous through the origin and defined on the real interval 
[— l/e,oo) is called the principal branch and is denoted 
Wq. The secondary branch, conventionally chosen to be 
the one defined on [— 1/e, 0] and denoted W-i, diverges 
at the origin (such that W-\ — > — 00). In our case, we 
have to choose the principal branch since for 6=1 the 
argument of W is positive, and for b = — 1 we must have 
< I. 

In the case of the running- mass model H35fl . the total 
number of c-folds can also be obtained explicitly. It reads 
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where the exponential integral function is defined by |24| 
Ei(x) = -J^dte-t/L Obviously, this expression is 
too complicated to be inverted. However, if, as done in 
Ref. |2l|, one notices that f c \/M pl -C 1 then one can just 
replace V in the expression giving the number of e-folds 
by M 4 . This leads to 



N~ - | In 



In^i 


-In 


ln^ 


) 






fa 





(41) 



and then, since the previous expression can be inverted, 
one obtains the classical field as a function of the number 
of e-folds explicitly, namely 



<p ci {N) = V > exp e civ lnrifi 



(42) 



Let us notice that the expression (14 1|) is in agreement 
with, for instance, Eq. (21) of Ref. |21| . 

Our next move is to find the numerical value of the 
parameters M and \i or c and ipo- This can be done from 
the measurement of the CMB anisotropy made by the 
Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite, namely from the formula 



Qrms-PS 

5C 2 



1 H 2 



J^2 



47T 607re* m 2 45e* m 4 



(43) 



where Q lms -Ps/T has been measured to be Q lms -Ps/T 
6 x 10~ 6 . We now discuss the four cases separately 



B. Large field models 

For a = we can eliminate the free mass parameter 
fi by simply rescaling the other parameter M. We thus 
obtain a monomial potential given by 



V{ip) = M 4 



(44) 



In this case, the slow-roll equation of motion leads to a 
solution which is completely explicit and reads 




47T 



-N . 



(45) 



The total number of e-folds during inflation is simply 
given by jV t = Air{tpi n /m pl ) 2 jn — n/A and can be very 
large if the initial energy density of the inflaton field is 
close to the Planck scale m 4 . The model remains under 
control only if the initial energy density is smaller than 
m pl and this imposes a constraint on the initial value of 
the field, namely ipi n /m pl < {m pl /M) 4 / n . 

For this kind of potential, the slow roll parameter 
e ~ m pi / {16n)(V /V) 2 becomes e = n 2 / (m pl / ip) 2 
and inflation stops when e = 1, i.e. when the slow- roll 
conditions are violated. The corresponding value of the 
field is (fend — n/(4y / 7r)m pl . The classical solution (|45|) 



allows us to calculate the field value <p* at Hubble cross- 
ing during inflation in terms of A* , the number of e-folds 
between the Hubble radius crossing and the end of infla- 
tion. We get (ip^/m pl ) 2 = n(n + 4JV»)/(167r) from which 
we deduce the corresponding value of the slow-roll pa- 
rameter of e* = n/(AN* + n). Finally, from the WMAP 
normalization, see Eq. (|43|l , we deduce the mass scale M 



M 



(45/2)n 



{AN* + n) n l 2 + x 



16tt 

n 



»/2 |Q2 



rms-PS 
J^2 



(46) 



This is the value of M that we use for the calculation 
of the quantum effects. From an observational point of 
view, all the models such that n > 5 are now excluded by 
the WMAP data, the quartic case being on the border 
line E3. 



C. Small field models 

In this subsection we discuss the WMAP normalization 
for the potential Ij34|l with a = 1, b = — 1 and n = 2. For 
such a model the slow roll parameter e reads 



Ann 2 [1 - {v/») 2 



(47) 



and, imposing e{<fi cn d) = 1 at the end of inflation, we can 
solve for cn d and obtain 



fond 

n 



16-KfJ, 2 



(48) 



Moreover, at Hubble radius crossing, the value of the field 
can be expressed exactly as 



-W 



^cndJ c -(y cnd / A1 )^m ri jV,/(27rp 2 ) 



(49) 



At this stage, it is interesting to introduce a second slow- 
roll parameter, £2, defined by e 2 = dlne/diV j2(|. Then, 
the spectral index can be written as n s = 1 — 2e — £2, 
where the slow-roll parameters are evaluated at Hubble 
radius crossing. Working out explicitly the above formu- 
las, one arrives at 



m 2 pi 1 + 



2nfi 2 



1 - {ip*/ny 



(50) 



Then, one can take the following route. If one chooses 
a value for the scale fx, then one can calculate ip cn d with 
Eq. |@SJ, ip* with Eq. and, finally, the spectral in- 
dex with the previous formula. Hence, instead of working 
with /1, one can express everything in terms of n s . Fi- 
nally, using Eq. (|43|l . one can determine the scale M in 
terms of ip„ and fi or, equivalently, in terms of the spec- 
tral index. In other words, we end up with an exact 
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relation M(n s ). It is represented in Fig. ^ (solid blue 
line). As one can notice on this figure, the curve blows 
up at n s = 1 — 6/(1 + 4A»). Let us try to understand 
this behavior in more details. If we expand the expres- 
sion giving the spectral index in terms of the parameter 
rn pl //i, one obtains n s = 1 — 6/(1 + 4A*) + m pl /[//(l + 
47V H ») 2 y / 7r] + O (m pl //i) from which we can obtain an ex- 
pression of m pl //i in terms of n s . In the same manner, 
one can expand M and the result reads (M/m pl ) 4 = 
45^ 2 /[(l+4Ag 3 / 2 (m pl / Ai )] (Q rms _p S /T) 2 +0[(m P1 / M ) 3 ]. 
Putting these two formulas together, one finally gets that 



M 



X ( 1 - n a - 



Qrms-PS 



T 



6 



1 + 47V* 



(51) 



up to terms of order O [(m pl //x) 3 ] , and we now under- 
stand the presence of the singularity. Although exact for 
potentials of the form (|34l) (with a = 1 and b = —1) , 
the behavior of M is this regime is not realistic for the 
following reason: on general grounds, it is clear that it 
is necessary to consider additional terms in the poten- 
tial l|34|) because, otherwise, tp = fi is not a minimum. 
If, during slow-roll inflation <p <C /i, then these terms are 
not important for the calculation of the perturbations. 
But, if fi 3> ti p1 , then one has that 4> C nd ~ A* and we 
expect the extra terms to play a role even during slow- 
roll inflation. In this regime, the shape of the potential 
that we used in order to obtain that M has a singularity 
is therefore not realistic. Moreover, it has been shown 
in Ref. |27| that the presence of these extra terms can 
strongly modify the energy scale of inflation. 

Let us now try to understand the behavior of the curve 
far from the singularity. This can be analyzed in the 
regime where </>*//i -C 1. In this situation, Eq. (|5U[) tells 
us that 



1 



(52) 



Moreover, since the Lambert function is small for small 
values of its argument, this also implies from Eq. i|49|) 
that [where we use Wq(x) ~ x] 



fend 



exp 



2>cnd 



2 N* ml 



2tt h 2 



(53) 



The previous derivation is consistent as long as the term 
proportional to A» in the argument of the exponential 
in i|53|) is large rendering the argument of the Lambert 
function small. This implies that A*m pi /27r/i 2 > 1 and 
provides the consistency constraint ( 1 — n s ) > 1 . This 
means that, in order for our approximation to hold, n s 
cannot be too close to 1 which is fine since this is precisely 
the regime that we are interested in [as we are trying to 
approximate M(n s ) far from the singularity]. Another 
way to see the same thing is to remark that, if n s were 



0.002 



0.0015 



e o.ooi 



0.0005 





\ 


\ 

\ 

\ 
















\ 




/ 












\ 





















0.1 



0.2 

l-n s 



0.3 



0.4 



FIG. 1: Characteristic mass scale M obtained from the 
CMB normalization as a function of the spectral index for the 
new inflationary scenario. The exact, numerically calculated, 
curve (solid blue line) blows up at 1 — n s = 6/(l + 4iV») while 
the approximated one (red dashed line) vanishes at the ori- 
gin. A different approximation, usually found in the literature 
(see for instance Refs. |2§l l2Sjp. is also shown for comparison 
(green dotted-dashed line). The difference between these two 
lines is approximatively a factor of 2 in M and is significant 
since this leads to a factor of 16 in the normalization factor 
M 4 of the potential and, thus, in the variance and the mean 
value of the fluctuations. 



too close to 1, then, from Eq. (J52J, it would follow that 
/i m m , in contradiction with the hypothesis that is 
exponentially damped, see Eq. (|53|l . 

One the other hand, as is apparent from Fig.^ we are 
not interested either in large values of 1 — n s which are, 
anyway, observationally excluded. Therefore, a Taylor 
expansion in this quantity is still valid in our case. In 
particular, using Eq. I|52l) . one can express (p cn d, hence 
fit and e*, in terms of 1 — n s only. Finally, since V* ~ M 4 , 
expanding everything in 1 — n s , we can write Eq. I|43|l as 



45 Q r 



^-^(l-nje- 1 -^ 1 -^), (54) 



T 2 



and we recover exactly the correct damping behavior ob- 
served in Fig. (red dashed line). 

Usually in the literature [28l |29| |. in the regime 
/x/to p1 <C 1, one expands Eq. H48J) to arrive at the expres- 
sion (ifend/n) ~ 2y / 7r/i/m pl . Then, following the same 
steps as before, we end up with 



M 



45 Q r 



2 

ms-PS 



-N,(l-n ) 



T 2 



(55) 



This expression has to be compared with Eq. JSJ 
The corresponding curve is represented in Fig. ^ (green 
dashed-dotted line) . For n s — 1 not too close to zero, the 
two predictions are in good agreement. 

Based on the above considerations, we arbitrarily 
choose to work with n s ~ 0.93 (compatible with the 
WMAP data) for which it follows that \i ~ 1.5m pl . In 
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this case, Fig. ^ indicates that 



M 



10" 



(56) 



This value corresponds to a regime which is compatible 
with the various approximations discussed above. We 
will adopt these values when plotting the quantum effects 
that we calculate in the following. 



D. Hybrid inflation 

This case is slightly different since hybrid inflation is 
in fact a two-fields model with the potential 20] 



V {y, V) = \™*<f + ~ - A 2 ) 2 + , 



(57) 



where ip is the inflaton and ip the waterfall field. A' and 
A are two coupling constants. The advantage of hybrid 
inflation is that inflation can be realized even if these 
coupling constants are of order one. The inflationary 
valley is given by ip = and, in this case, the potential 
reduces to a potential similar to the one given in Eq. 134|) 
with a = 1 and 6=1, provided we have 



M = 



A ^/ 4 A 
V2 



A' A 2 



(58) 



In the hybrid scenario there are a priori two mecha- 
nisms for ending inflation. Either inflation stops by in- 
stability when the inflaton reaches a value 



(59) 



where the mass in the direction perpendicular to the in- 
flationary valley becomes negative or the slow-roll condi- 
tions are violated which happen for 



\ 




(60) 



Let us notice that (f> e does not exist if fi > m p J '(4^/tv). 
Therefore, the final value of the inflaton, cn d, is the 
maximum of <^ cr ; and <p e . To decide which mechanism 
is realized in practice requires the knowledge of the pa- 
rameters of the model. Of course, the parameters of the 
models must satisfy the CMB normalization [2(J. The 
value of the field at horizon exit is given by an equation 
very similar to Eq. (|49|l . namely 



\ 




exp 



- 1 



(61) 



Then, one can repeat the algorithm described in the 
previous subsection to implement the CMB normaliza- 
tion. Here, we assume that the coupling constants are 
A = A' = 1 and use the following values 

-8™ 



A = 10" 



in 



0.7 x 10" 



(62) 



This implies fi ~ 1.01to p1 and n s ~ 1.08, i.e. a blue 
spectrum as expected for hybrid inflation. Let us remark 
that in the above case, inflation ends by instability and is 
vacuum dominated. This is the regime of interest for us 
because it would be pointless to calculate the quantum 
effects in the regime where the field dominates since this 
case reduces to the case of chaotic inflation. 



E. Running-mass Inflation 

Running mass inflation can be realized in four differ- 
ent ways [21| that we now very briefly describe. In the 
following, we refer to these four different models as RM1 
to RM4. From the expression of the potential (|35|l . it is 
easy to see that (po is an extremum of V(<p). This is a 
maximum if c > and a minimum if c < 0. According 
to the classification of Ref. 21], "model 1" (RM1) cor- 
responds to the case where c > and p c \ < pa- In this 
case, ip c \ decreases during inflation. "Model 2" (RM2) 
also corresponds to c > but, now, with p c \ > tpo and 
ip c \ increases during inflation. "Model 3" (RM3) refers 
to the situation where c < 0, for which tpo is a minimum, 
and p c \ < tpo all the time. In this case, ip c \ increases dur- 
ing inflation. Finally, "model 4" (RM4) has c < and 
Pel > <po decreases during inflation. 

The values of the free parameters c and ipq are con- 
strained by the CMB measurements. In Ref. |2l|. these 
constraints have been studied for the four models evoked 
above in the parameter space (c, a) , where a is defined 
by a = — cln ((/?end/Vo)i Pcnd being the value at which 
inflation stops. Using a allows us to express the spectral 
index and the running of the spectral index as 



1 



-2c+ 2crc 



-cN. 



= -2c- 



1 dn„ 



c d In k ' 



(63) 



where iV* has already been defined before. 

Let us now study how inflation ends in those models. 
A priori, the end of inflation is found from the condi- 
tion e = 1. However, it is easy to see that this cannot be 
achieved for model 3 [for V((po) ^ 0] and model 4. In this 
case, another mechanism must be advocated, presumably 
of the hybrid type. For simplicity, in the following, we 
focus on models 1 and 2 only. In the case of model 1, 
one can also show that the condition e = 1 cannot be 
satisfied because e — > when ip — > and is bounded by 
c 2 ip^l (l&neM p ) <C 1 in the interval [0, p>o\. However, the 
slow-roll parameter r\ = M 2 V" /V (this slow-roll param- 
eters differs from £2 introduced above) increases when 
ip — > and inflation stops when X] = 1. In the case of 
model 2, the condition e = 1 is a priori possible but we 
will sec that, in practice, the condition r\ = 1 occurs ear- 
lier and, therefore, controls the end of inflation as for 
model 1. 

As an illustration of model 1, we work with the values 
c = 0.06 and <po/M pl = 10~ 6 . In this case, the end of 
inflation occurs at ip en d/M pl ~ 2 x 10~ 14 . The param- 
eter a is given by a ~ 1.064. One can check in Fig. 3 
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of Ref. [21j (left panel) that these values are compati- 
ble with the CMB constraints and, moreover, that they 
correspond to a situation where the Planck satellite will 
measure the running of the spectral index. 

For model 2, the situation is slightly more complicated. 
If we use the same parameters c and ipo, then inflation 
ends at ip en d/M pl ~ 0.6 which is fine since this value 
is smaller than one. However, this leads to a ~ —0.79 
which, according to Fig. 4 (left panel) of Ref. [2lJ is not 
acceptable. Clearly, a way to cure the previous problem 
is to increase the value of Pq. Therefore, for instance, 
we could try ip /M pl =0.1 and c = 0.06. In this case 
<Pend/M P1 ~ 2.95 and a ~ —0.2. The values of c and 
a are now compatible with the CMB constraints but, 
of course, the value of p C nd > M pl is problematic. For 
instance, the simple expression of the number of e-folds 
given by Eq. (|41|l would not be valid in this case because 
we have assumed p c \/M pl <C 1. Moreover, we should 
also include higher orders terms in the potential. For 
these reasons, and since we only want to illustrate how 
quantum effects behave in the running mass scenario with 
simple formulas, we will continue to assume that c = 0.06 
and ipo/M pl = 10^ 6 even for model 2. 

Finally, once c and po have been chosen, the COBE 
normalization fixes the scale M through the relation 



M 4 



45 <9 rms -ps v 2 <fl ( cr _ cAr 

t, 7T- exp e *- 

m- 4 T 2 Ml V c 



2dV„) . (64) 



In the case of model 1 this gives M ~ 6.4 x 10 



while 



m pl , where in both 



for model 2 we have M ~ 1.2 x 10~ 6 
cases we have used iV» = 50. 

Having chosen the parameters, one must also specify 
the initial conditions. The total number of e-folds is given 
by 



1 



In 



In 



Pin 
Po 



111 



Pcnd 
P0 



(65) 



Regardless of the sign of c, one can check that N T > 
0. For model 1, one takes p- m = (l — 10~ 5 ) i^o which 
gives N T ~ 240, while for model 2 one chooses ipi n = 
(l + 10~ 3 ) <po which implies that N T ~ 165. 



RESULTS AND DISCUSSION 



with 



J„ = 4 r n dxx z{1 - n) (a + bx n f 



= Pn 



fin 



-P» 



fc]_ 



the function P n being given by 

„ , , 4 /, 4a x~ 3n 6abx- 2n 
P n (x) = x 1 6 + 



4- 3n 



n-2 



12ax~ n 

71-4 



(67) 
(68) 

(69) 



This expression is valid as long as n ^ 2,4 or 4/3. These 
three particular cases must be treated separately and the 
corresponding function P reads 



Pi{x) 



-2^r + 6ax z 
ar 

ab 



bx + 12abln:E . 



- 3— + bx 4 + 12a In a; , 

2x s 

P 4/3 (x) = 9abx i/3 + ^-ax s/3 + bx A + 4a\nx . 



(70) 
(71) 
(72) 



It is worth noticing that setting a = and b — 1 we 
correctly recover the result obtained for a simple power 
law potential V(tp) oc <p n , namely 



(difij) _ 4 M 4 

u 2 3n m 4 
^ pi 



gel 



n-2 



fin 
/' 



Pel 

A* 



(73) 



This result was also derived in Ref. [13 ■ 

The next step consists in calculating the correction to 
the mean value. Using Eq. (|15fl . this is immediately done 
once the variance is known and the result is 

M-n M) 1 AbMi 

~^ n[ ' + 3nm 4 



, (Pd/vT- 1 



+ b((p c i/fj,y 



fin 
fl 



FL 



Pel 

A 4 



(74) 



where the two functions Q n and i?„ are defined by the 
following expressions 



Qn(x) 



2a(n- l)+b(n-2)x n 
4x(a + bx n ) 



(75) 



Having fixed the values of the parameters for the four 
models under consideration, we can now compute the 
quantum effects. Let us first apply our formalism to 
the potential Ij34(l and study the behavior of the fluctua- 
tions in the different inflationary scenarios. In this case, 
Eq. (|T3|) can be integrated exactly and straightforward 
calculations lead to 



46 M 4 



3n to 4 a + b((p c i/fj,) n \ fi 



Pd 



2(n-l) 



(66) 



and 



R n {x) = x 2[1 - n) {a + bx n f /2 



(76) 



Let us notice that, for a = 0, these functions are simple 
power-laws and that, in this case, we recover the formulas 
found in Ref. [13j for generic monomial potentials. 

The last step consists in computing the volume effects. 
For this purpose, we have to compute the integral in 
Eq. for the potential This leads to 
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(<p) v -(<p) _ 32ttMV 



'Pel 



^a + b(<p cX /n) n 



(77) 



where S n is the function defined by the expression 

w = / dx (a + ( b :;f 2 (78) 

7/2 ^/4-3n 7 4-2n 6 \ 
= <* ' i — ^2^1 ; — x n ), 

4 — An \ n I n a J 

for 0. In this expression 2^1(0, 6; c; x) is the hyper- 
geometric function taking real values for values of the 
argument x less than one 0. This is indeed the case 
since, for 6=1, the argument is negative and for b = — 1 
we always have (ip c \/fj,) n < 1. If a = (in the large field 
models case), then the argument is ill-defined. However, 
using the relation for the hypergeometric function of ar- 
gument 1/x [2^| (that can be safely applied since, in this 
case, we have b = 1) we obtain a well-defined expression 
that yields 

T 4+n/2 

S n (x)=2-—, (79) 
8 + n 

up to an unimportant additive constant. 

Once the variance and the mean value are known, 
calculating the probability distribution of the field is 
straightforward, see Eq. I)25|l . Results for large field and 
small field models are shown in Fig. [2J For LF, we take 
an initial condition ip- in such that the potential energy is 
close to the Planck scale m* . In this case, P c (or rather 
its maximum (p) c ) slowly rolls down the potential while 
remaining "behind" the classical solution. At the same 
time and very quickly after the initial time, the variance 
significantly increases, i.e. P c strongly spreads, and, as 
a consequence, the tail of P c penetrates into the region 
where <p < 0. On the contrary, P v , while also spreading a 
lot, immediately inverts its motion, starts rolling up the 
potential and penetrates into the trans-Planckian regime 
where V . These behaviors confirm the impor- 

tance of quantum fluctuations in the LF models but this 
also rises the question of the reliability of the solutions 
obtained before. Indeed, in the case of P c , for ip < ip en d (a 
few m pl 's) the Hubble crossing of the modes stops and the 
noise no longer exists. Moreover, the slow-roll approxi- 
mation that we have explicitly used in our derivation of 
the Langevin equation is no longer valid. With regards 
to P v , in the regime where V > rnc , the whole frame- 
work of quantum field theory itself breaks down. These 
features are also very sensitive to the initial conditions. 
For instance, if we have initially V(ipi n ) < m^J 10, then 
the stochastic deviations from the classical trajectory are 
already dramatically reduced. 

In the case of SF models similar conclusions hold. If 
Pi n is very small (close the maximum of the potential), 



then P c significantly spreads around its mean value while 
rolling down the potential. The only difference is that, 
now, the peak of the distribution stays ahead of ip c \. As 
before, the tail of the single-point distribution goes to the 
region where the slow-roll approximation breaks down, 
P > fend and p > u. In this case, the previous con- 
siderations on the reliability of the solution still apply. 
Let us now study the behavior of P v . As can be seen in 
Fig. [21 at the beginning of inflation, the behaviors of P c 
and P v are similar but, after a few e- foldings, P v reverses 
its motion and starts moving back towards the maximum 
located at ip — 0. After enough time, most of P v is on the 
other side of the potential. Diffusion is less crucial than 
in the LF models because the energy scales involved are 
smaller but, nonetheless, we see that such effects are very 
important. On the other hand, the behavior of P v is less 
problematic since the Langevin equation is fully trustable 
even if the field is negative (provided, of course, it is 
still small in comparison to fj,). The only danger comes 
from a possible breakdown of the perturbative expansion 
as the field goes to a region where the stochastic mean 
value is very far from its classical counterpart. In this 
case, the quantum effects should in fact be computed as 
perturbations of a classical solution living in the region 
ip < 0. Then, the volume effects will act the same way 
as before, that is to say they will push the correspond- 
ing distribution back to the origin, towards the region 
p > 0. Therefore, it seems reasonable to postulate that 
a situation of dynamical equilibrium could set up with a 
stationary solution for P v concentrated around tp = 0. 

Finally, let us notice that, as before, the stochastic be- 
havior of the field is strongly dependent on the initial 
conditions. If ip is far from the maximum, then the drift 
due to the classical term dominates on the noise-induced 
diffusion, the effect of which is therefore no longer im- 
portant. The dependence on the initial conditions can 
be also be checked with the help of Eq. (|66|l in the limit 
Pel "C /i. In this regime we have ^2(2;) — — 2/a; 2 and, 
since ip- m < p c \, the relative amplitude of the fluctua- 
tions can be written as 

vm<ji(iL) 2 jL„ 10 -i, (80) 

Pel V 3 \ m pi/ fin 

for M ~ 10~ 3 m pl and p- ln ~ 10~ 5 /i. This shows that, 
in this case, the stochastic fluctuations can represent up 
to 10% of the field value (and they increase in the sub- 
sequent evolution since, in new inflation, different clas- 
sical solutions diverge). However, taking another initial 
condition, for instance <pi n /M — 10~ 3 , would reduce the 
fluctuations to only 0.1% of the classical field. Therefore, 
we see that there is indeed a strong dependence on the 
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initial conditions. 

Let us now consider hybrid inflation in the vacuum 
dominated regime. In this case, the corresponding distri- 
butions are so peaked that they cannot be distinguished 
from a Dirac distribution. This is the reason why we 
have chosen not to include any figure for the hybrid case. 
The reason for this behavior can be easily understood if 
we come back to Eq. (|66|) . In this case, ipi n > p c \ > p cr i 
and, therefore, one has 

ip c \ V 3 \ TO Pi/ 

the upper limit of the above expression being obtained 
when evaluated for ip = p cr [. Using Eqs. (|58(l and (|59|l . 

one arrives at 

^><-*i-~10- (82) 

for the values of A and m chosen before. Therefore, in 
hybrid vacuum dominated models, the fluctuations are 



less than 0.1% of the classical field until the very end of 
inflation. In addition, if one uses Eq. (|74|l in the same 
limit, one sees that the two terms of that formula ex- 
actly cancels out and (8^2) ~ 0. This means that, with 
a very good accuracy, the distributions are peaked over 
the classical value of the field. This explains the behavior 
of the distribution functions described before. Moreover, 
this conclusion is rather independent of the initial condi- 
tion, provided that the latter does not violate the vacuum 
dominated regime. 

Let us now investigate the quantum effects for the two 
different models of running mass inflation. The variance 
and the mean value can be calculated exactly since the 
corresponding integrals are feasible. This results in quite 
complicated formulas which are not especially illuminat- 
ing. However, since the vacuum expectation value of the 
inflaton is always small with respect to the Planck mass, 
it is therefore a good approximation to replace V /V by 
V'/M 4 as we have already done before, see Ref. [21| . 
Then, the variance reads 



Ml 12n 2 c\M p J Uo/ Wo 



1o 

Pin 



1 
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<Po\ 1 
p ci J ln(<£ c i/Vo) 
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-2Ei 
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^From this result, using Eq. (|15|l . we deduce the expression of the mean value, namely 

1 + In (ifd/ipo) , c p <£ci \ (8<pT) ' ' ' ' " ' 



(8p>2 



1 



fcci/MpOMVciM)) 2M P , ^ J Ml 
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48tt 2 c I M. 



M 



'Pel , / <Pd 

M FI n \p Q 



<Pcl ) In 2 (pci/tpo) 



(83) 



(84) 



Finally, let us compute the volume effects given by 
Eq. H33|) ■ Following the approximation already used 
above, one could take a factor H out of the integral. 
Then, the remaining integral, the kernel of which is now 
(H/H 1 ) 3 , can be performed explicitly. However, in this 
case, the result exactly cancel the second term in Eq. Q33JI 
and we would obtain (tp) v — (ip) ~ 0. Let us notice that 



this is is not a specific feature of the running mass po- 
tential but a general consequence of the approximation 
used. Therefore, in order to compute the volume effects, 
one should evaluate the original integral. Unfortunately, 
an exact integration is not possible in this case. A pos- 
sible solution is then to expand the integrand in powers 
of c(ipo/M pl ) 2 . At first order, this gives 



3I T J 



7M M 3 
16tt 2 c M p 3 M ; 
1 



'Pel , Vcl 

In — 

<Po 



1 - 



2 V 2 



A. 

Ml 



-1/2 



1 



h\(p- m /p a ) ln(y c i/Vo) 41n 2 (< ( £>in/< ( fo) 41n 2 (y) c i/^o). 



(85) 



Alternatively, a numerical integration of Eq. (|33|) can al- ways be done and, in the following, we present results 
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FIG. 2: Evolution of the single-point (solid blue line) and volume-weighted (dashed red line) probability distributions, P c (<p) 
and P v (</>), for the large field model V oc <p 2 (LF) and the small field model V cx 1 — (ip/fi) 2 (SF). The initial values are 
(pin = 7.6 x 10 5 m pl (corresponding to Vi n = m p ,/2) for LF and (pi n /ft = 1CF 5 for SF. The initial shape of the two probability 
density functions is always chosen to be 8(tp — <pi u ). The vertical dotted black lines represent the location of the classical field. 
Three successive snapshots of the distributions (numbered 1, 2 and 3) are shown on the left panels while the evolution of (p) c 
and (<p) v is displayed on the right panels. The classical field ip c \ evolves from the right to the left in LF and from the left to 
the right in SF. In both cases, P c {<p) rolls down the potential, spreads significantly around its mean value and penetrates into 
a classically forbidden region (ip < for LF and ip > fi for SF). The quantity {tp) c stays "behind" the classical value in LF 
but is "ahead" in SF. On the other hand, P v (<p) (not shown in the LF left panel) reverses its motion and climbs towards the 
trans-Planckian region (in LF) or towards the maximum of the potential at ip — in SF. 



obtained with this method. 

Results for RM1 and RM2 are displayed in Fig. EI Let 
us start with the RM1 model. In this case, P c {<p) fol- 
lows very closely the classical solution, the peak of the 
distribution being slightly "behind" <p c \ (as for LF). In- 
terestingly enough, if, at the beginning of inflation, the 
distribution P c {<p) starts spreading around its mean value 
as it was the case for the LF and SF models then, after 
some e-foldings, the variance reaches a maximum and 
then starts decreasing, i.e. P c (<p) becomes more and more 
peaked over the classical solution. This is a consequence 
of the fact that the classical dynamics (which dominates 
at late times) tends to attract different solutions contrar- 
ily to the SF case where they are instead pulled apart. 
The behavior of P v (<p>) is even more interesting. The 
volume weighted distribution moves backwards, beyond 
the maximum of the potential, into the region classically 
corresponding to RM2. Then, it reverses its motion and 
comes back into the region corresponding to RM1. This 



last behavior is probably not trustable because, as long as 
P v (<p) penetrates the region RM2, the calculations per- 
formed before should be modified to take this situation 
into account. Let us also notice that this behavior is 
not exactly similar to the one observed for the SF model 
where the distribution goes to the region <p < 0, see the 
discussion above. Indeed, in this last case, whatever the 
sign of the field, the model is the same, in particular the 
value of M remains unchanged. On the contrary, RM1 
and RM2 are really two different models with two differ- 
ent energy scales. Let us now study RM2. As for RM1, 
the peak of P c (<fi) follows the classical solution with the 
difference, however, that it is slightly ahead <p c \ (as for 
SF) and that the spreading of the distribution continu- 
ously increases. At the beginning of inflation, P v ((p) also 
follows <p c \ but, then, it changes its motion and moves 
back towards the region corresponding to RM1. 

One of our main conclusion concerning running mass 
inflation is that, if RM1 and RM2 are classically two dif- 
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FIG. 3: Evolution of the single-point (solid blue line) and volume-weighted (dashed red line) probability distributions for the 
two different running mass inflation models RMl and RM2. The initial probability density function is S(<p ~ ipi n ) for the two 
models. The initial values are chosen to be (fi a /<po = 1 — 1-5 x 1CP 5 for RMl (<p < ipo, inflation proceeding from the right to 
the left) and v?in/Vo = 1 + 10 -3 for RM2 (p > 950, inflation proceeding from the left to the right). The vertical lines (dotted 
black lines) represent the location of the classical field. On the two left panels, three snapshots of P c (<p) and P v (p) (numbered 
from 1 to 3) are shown at three different times (respectively corresponding for RMl to p/po = 0.9, 0.5 and 0.01) together with 
the corresponding values of the classical field. In the right panels, the evolution of (ip) c and {tp) is followed until the end of 
inflation. The physical interpretation of these results is discussed in the text. 



ferent models, from a statistical point of view, it may 
then be impossible to distinguish between them. Indeed, 
regardless of which model we decide to start with, if the 
initial conditions are sufficiently close to (fo, then there 
is a significant probability of diffusing on the other side 
of the potential. Moreover, when we start with RMl 
(RM2), volume effects tends to push the evolution to- 
wards RM2 (RMl). Therefore, as it was the case for 
SF, it seems reasonable to postulate that the system will 
settle down at the boundary between the two models, a 
situation described by a stationary distribution concen- 
trated around ipo. Such a situation typically leads to the 
sclf-rcproducing regime. 

Finally, another important feature of the running mass 
potential is that quantum effects can strongly modify the 
classical evolution even if the energy scale involved is far 
below the Planck energy (this is also the case for new 
inflation). This reinforces the fact that the connection 
between the importance of the quantum effects on one 
hand and the fact that V is close to the Planck energy 
on the other hand is very specific to the large field mono- 



mial potentials used in chaotic inflation. In fact, it is 
clear that quantum fluctuations are important whenever 
the classical contribution to the motion of the field is 
suppressed. In chaotic inflation, this happens close the 
Planck scale because of the large friction term but this 
also happens near the maximum of the potential (and far 
from the Planck scale) in new and running mass inflation 
models because of the smallness of V'(ip). 

We end this section by a slightly different discus- 
sion. In Ref. eternal inflation is studied using the 
Langevin equation for the model V(<p) cx ip 4 . The main 
argument presented for considering this simple model is 
that analytical progress may be made in this case. Us- 
ing the method presented in our article, one can in fact 
consider a much larger variety of models and one is not 
restricted to the simple chaotic quartic model. In order 
to illustrate this claim we consider the calculation, pre- 
sented in Ref. 11], of the value of the inflaton at which 
the self-reproducing behavior becomes possible, not only 
for the potential <p 4 but for the general case ip n (in fact, 
one could even reproduce this calculation for new infla- 
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tion, run ning mass inflation and so on). 

InRef. [11|, the main idea is to evaluate the mean value 
of the number of e-folds is given by 



(N) = 



dt (H) = 



<p(T) 



(86) 



This expression can easily be evaluated for any potential 
using our perturbative treatment. For large fields models, 
straightforward calculations lead to 



(AO = N* 



1 



4(3n + 4)M 4 
3n(n + 2)m 4 



(87) 



The breakdown of this expansion signals the beginning 
of the self-reproducing regime. Using the fact that A^ 1 ~ 
47r(yi n /m pl ) 2 /n, one easily sees that this happens when 
the initial value of the field is 

jftn „ , (88) 

m pl 

where the coupling constant A„ is defined such that 
M 4 /m™ = A„/to™~ 4 . For n = 4 one recovers the condi- 
tion found in Ref. Our method allows us to obtain 
this condition for any potential and there is no need to 
assume a quartic potential to perform this calculation ex- 
plicitly. Incidentally, the above condition was also found 
previously for instance in Ref. Q [see Eq. (3.31)]. 



equation rather than in its solution. In particular, we 
were able to provide a general second order expression 
for the probability distribution of the field with or with- 
out the volume effects taken into account. Our expres- 
sion only requires the calculation of one quadrature (two 
if the volume effects are considered) which, for most of 
the inflationary scenarios, can be performed explicitly. 
Second, we have applied this method to various models 
of inflation. This has allowed us to compute the quan- 
tum effects with backreaction for chaotic, new, hybrid 
and running mass inflation. To our knowledge, in the 
case of the last three models, this is the first time that 
such a calculation is done (the backreaction being taken 
into account). Third, we have discussed the impact of 
the stochastic effects on these inflationary scenarios. For 
instance, in the case of running mass inflation, it was 
shown that the quantum effects blur the distinction be- 
tween the various running mass inflationary models and 
that the self-reproducing regime is likely to be important. 

Finally, an important advantage of the method pre- 
sented here is that its accuracy and domain of validity 
can be evaluated in details. This will be the subject of a 
forthcoming paper [l5| . 
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APPENDIX A: GENERAL DEFINITIONS 



In this short Appendix, we give the precise definitions of the quantities that we have been using in order to establish 
the expression of the probability density function in Sec. IIIII C. If we define a "two-step" function l d{t- m < r < t) = 
i?(t — t)i9(t — t in ) (whose value is 1 if the inequality is true and otherwise), the expression of the "vector" J is 



Jt.tjr) _*(«,, <r<t) ^ 



(Al) 



while the "matrix" A can be written as 



(ci,ct 2 ) 



H' cl (t) 

47T 



3tf(t in < ai < i)V ' H c i{ai)J ait t ln {a 2 ) 



^r,t in (o'l) J r , 4in ((T2) 



(A2) 
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Finally, we also provide the definitions of quantities that have been using in Sec. IHII D where the volume effects have 
been estimated. In particular, the "vector" / can be expressed as 



Jm„» =3jf dTH' cl (T)J T . tia (a), 
and the "matrix" B is defined by the following formula 



Bt itin (eri,cr 2 ) = jf dr 



Kir) 



(A3) 



(A4) 
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